Moving boundary approximation for curved streamer ionization fronts: Solvability 

analysis 
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The minimal density model for negative streamer ionization fronts is investigated. An earlier 
moving boundary approximation for this model consisted of a "kinetic undercooling" type boundary 
condition in a Laplacian growth problem of Hele-Shaw type. Here we derive a curvature correction 
to the moving boundary approximation that resembles surface tension. The calculation is based on 
solvability analysis with unconventional features, namely, there are three relevant zero modes of the 
adjoint operator, one of them diverging; furthermore, the inner/outer matching ahead of the front 
has to be performed on a line rather than on an extended region; and the whole calculation can be 
performed analytically. The analysis reveals a relation between the fields ahead and behind a slowly 
evolving curved front, the curvature and the generated conductivity. This relation forces us to give 
up the ideal conductivity approximation, and we suggest to replace it by a constant conductivity 
approximation. This implies that the electric potential in the streamer interior is no longer constant 
but solves a Laplace equation; this leads to a Muskat-type problem. 



I. INTRODUCTION 

Streamers are growing plasma channels that in suf- 
ficiently strong electrostatic field expand into ionizable 
matter like air or other gases, liquids or solids. The 
field accelerates the electrons sufiiciently to ionize gas 
molecules upon collisions. The self- focussing of the field 
at the tip of the propagating streamer strongly sup- 
ports this process. The richness of spatio-temporal struc- 
tures that form as streamers undergo successive branch- 
ing events during their propagation make them an impor- 
tant example of nonequilibrium pattern formation. They 
are of fundamental interest not only for this reason, but 
also because they determine the initial stages of electric 
breakdown in sparks and lightning, in technological ap- 
plications as well in natural processes [H, [3, S [j| • In this 
paper we focus on negative (anode-directed) streamers 
in simple gases like pure nitrogen or argon. Their dy- 
namics can be described by a minimal model of nonlinear 
reaction-advection-diffusion equations, as described, e.g., 
in [a, Ig] . The model describes the evolution of electron 
and ion densities and their nonlinear coupling to the local 
electric field. 

Many simulations 0, i, i, [13, [H El, [11 H since 
the early 1980'ies have shown that the ionized interior 
of the streamer finger is preceded and surrounded by a 
narrow curved space charge layer that enhances the elec- 
tric field in the non-ionized region ahead of the front and 
screens it in the ionized interior. In particular, in strong 
background fields after a sufficiently long evolution, the 
width of the ionization front can be much smaller than 
its radius of curvature [ij, [ij, [l^|. This separation of 
scales enables one to consider the front as an infinitcsi- 
mally thin, sharp moving interface F. Many inhomoge- 
neous systems involve domains of well defined phases sep- 
arated by thin interfaces. These include non-equilibrium 
systems like solidification occurring by dendritic growth 
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FIG. 1: On the left: generic solution (net charge density) of 
the minimal PDE model with curvature « and width l^ (see 
Eq. (|30|l 'l of the ionization front. On the right: the moving 
boundary approximation (MBA) with the ionized region, fl^ , 
the non-ionized region, f2+, and the sharp interface, F. 



[H, [13, [S [H, [13, [U, [H, bacterial growth 0, [H or 
many diffusive-reactive systems [2^. In Fig.[Tl we present 
a generic solution for the net charge density of the mini- 
mal PDE model showing the separation of scale together 
with the resulting moving boundary approximation. The 
original nonlinear dynamics is then replaced by a set of 
linear field equations (typically Laplace) on both sides of 
F, r2+ and f2^, with appropriate boundary conditions at 
the interface, r+ and F" respectively, and further away 
from it [2l|, [2y, [23, ^] (and references therein); the non- 
linearity enters through the motion of the boundary. The 
interface dynamics is then typically related to gradients 
of the Laplacian fields at its vicinity, F"*" and F~. 



In the context of streamer dynamics, the concept of 
an interfacial approximation was proposed aheady in the 
early 1970'ies by Lozansky and Firsov before any numer- 
ical simulations [29]. At the time, their model was in 
competition with other streamer models and concepts. 
Lozansky and Firsov proposed to consider the streamer 
interior, fi^, as ideally conducting, i.e., the electric po- 
tential (f) to be constant in the interior. The exterior, 
57+ , is non-ionized and therefore does not contain space 
charges, solving in electrostatic approximation (that is 
justified in ^) 



V 



in ^+. 



(1) 



The interface was assumed to move with the local elec- 
tron drift velocity 



v = V0^ 



(2) 



From here, superscripts attached to fields, potential 
and densities indicate their limit value as they approach 
the interface from O"*" and 17" , respectively. In particular, 
we denote (jy^ = 4>\y+ and (j)~ = 4>\y- ■ However, the 
solutions of this interfacial model were hardly explored. 
Only 30 years later, the interfacial model was taken up 
again in the streamer context |30i] , and Meulenbroek et al. 
showed that the model actually allowed for spontaneous 
streamer branching. However, the model was also found 
to be mathematically ill posed; in the context of fluid 
dynamics, this is explain for example in Ref. [3l| and 
references therein. To resolve this problem, the boundary 
condition 0+ — 0^ = was replaced by the regularizing 
boundary condition 



Qo(n • V0^ 



|n-V0+|>l , 

— > n ■ 



V(/>^ 



(3) 



which was proposed in [32l] and derived in planar front 
approximation in [33] ■ The boundary condition accounts 
for the finite width of the charged layer (the ionization 
front) that leads to a finite variation of the electric poten- 
tial across the front. The boundary condition in the limit 
of large electric fields actually turns out to be identical 
to the "kinetic undercooling" boundary condition that 
was applied to crystal growth under certain conditions 

[Mil. 

The solutions of the interfacial model with this bound- 
ary condition ([3]) implemented were studied in [32l |33J, 
l36l . l37j , and the results of these papers show that the in- 
terfacial model employed there is the simplest one that 
regularizes the motion - therefore we call this model the 
minimal regularized model. However, it also was shown 
recently that this analysis applies only to boundaries 
that at time t — Q are sufRciently many times differ- 
entiable [H. 

In this paper we compute curvature corrections to 
the boundary condition ([3]). A systematic expansion of 
slightly curved fronts about planar fronts for streamers 
was first suggested in Q , and Q contains the analysis of 
planar fronts as first step in this research program. In 



analogy with other weakly curved fronts, a perturbative 
expansion of the curved front about the planar front with 
subsequent solvability analysis was the method at hand. 
Solvability analysis means that the inner front structure 
is integrated out and replaced by conditions for inter- 
facial motion that are matched to the dynamics on the 
outer scale. 

However, such an analysis could not be performed on 
the streamer model of [a, [g| as the fronts are pulled; for 
pulled fronts, we refer to [39l,i4(j|, and for the impossibility 
of a solvability analysis for pulled fronts to 41]. The 
formal reason for the non-applicability of the method are 
diverging integrals in the leading edge of the ionization 
front; the physical reason is the algebraically slow front 
relaxation of the leading edge that actually stretches out 
through the whole non-ionized region and does not make 
part of the inner front region that is to be integrated out. 
The pulled nature of the front also requires special care 
with grid refinement in numerical solutions [14| . 

However, the leading edge that pulls the front, is diffu- 
sive, and it is a physically and mathematically meaning- 
ful approximation to neglect electron diffusion in strong 
fields, where electron motion due to drift dominates over 
the diffusive motion [i^, |43] . We have employed this ap- 
proximation in the interfacial model and we have checked 
that its results approximate simulations with electron dif- 
fusion well [33|, |4^ . We here employ the same approxima- 
tion which now makes it possible to perform a solvability 
analysis. 

In fact, neglecting electron diffusion and assuming that 
the state ahead of the front is completely non-ionized, 
converts the inner-outer matching problem from one non- 
typical situation to the other. With electron diffusion, 
the inner front region was not separated from the outer 
one. Without electron diffusion, however, the inner re- 
gion finishes precisely where the electron density discon- 
tinuously jumps to 0, and the matching between inner 
front region and the non-ionized outer region takes place 
precisely on this line rather than on an extended spatial 
region. 

The "moving boundary condition" (MBC) for a 
slightly curved front dynamics now can be systemati- 
cally derived from the original nonlinear field equations: 
A perturbation of a planar front is assumed whose cur- 
vature in the direction transverse to the front motion 
is much smaller than the front width, and solvability 
analysis is used to connect the perturbed values of the 
fields ahead and behind the curved front. Recently, such 
approach has been successfully applied to derive MBC 
equations for the dynamics of discontinuity curves which 
appear in nonlinear diffusion problems |24| . 

The solvability analysis poses technical challenges. We 
will show that the electrostatic field behind a slowly 
evolving curved front does not vanish but is rather con- 
stant and proportional to the front curvature. Conse- 
quently, the electric potential diverges linearly behind 
the front. Another divergence appears in one of the zero 
modes of the solvability analysis. This divergence is how- 



ever necessary to cancel the one of the electric potential 
behind the front in order to get a meaningful boundary 
condition for the jump in the electric potential across 
the interface T. To our knowledge, this nontrivial com- 
putational aspect is new, and we are not aware of similar 
types of solvability analysis in implementations of sharp 
interface approach to front dynamics in other physical 
systems. 

Furthermore, since the potential is no longer constant 
behind the front, the ideal conductivity approximation 
((/) = in f2^) must be relaxed and is replaced by 
a constant conductivity approximation for the ionized 
streamer interior (V^0 = in n~). 

Explicitly, our front dynamics consists of the following 
equations: 



with the MBC 



A(j) = in n+ 
A(j) ^ in n^ 



g2(n-V0+)K 

go(n-V0+) + Qi(n-V0+)K 

n- V(/)+. 



(4) 
(5) 



(6) 
(7) 
(8) 



where the coefficients Qi depend on the electrostatic field 
ahead the front, and are given by analytic formulas de- 
rived from the planar front solution, see Eqs. (|129p - (|13ip . 
The paper is organized as follows: In Sec.[lTlwe summa- 
rize previous results on the minimal model and its planar 
front solutions essential for our analysis. In Sec. lIII Al we 
start our formal derivation by assuming a slightly curved 
front with e = Ki~ <C 1, where £~ and k are the width 
and the curvature of the front respectively. We construct 
a perturbative expansion of the solution to the minimal 
model dynamics in powers of e around the zeroth order 
planar front solution. We focus on the leading order be- 
havior in e and describe it as a solution to a set of four in- 
homogeneous linear ODEs. In Sec. lIIIBl we explain and 
develop the solvability analysis formalism: Using zero 
modes of the adjoint linear differential operator we de- 
rive relations between deviations of the electrostatic field 
and ion/electron densities ahead and behind the curved 
front, r+ and T~ , from their values at the planar front so- 
lution. These relations turn out to be the required MBC 
equations for the curved fronts dynamics. The required 
zero modes are calculated analytically in Sec. IIII CI In 
Sec. IIIIDl we conclude our inner analysis and derive the 
MBC equations. The final form of the MBC equations 
(|124p - (|128p . obtained by matching inner and outer re- 
gions, is derived in Sec. IIVI where we give also the large 
field limit of the interfacial model together with an ele- 
mentary algorithm to solve it. In Sec. |Vl we discuss the 
constant conductivity approximation for the ionized re- 
gion. In Sec. IVII we conclude by highlighting the major 
results of this paper and point to future directions. 



II. COLLECTION OF SOME PREVIOUS 
RESULTS 

In this section we summarize previous results required 
for the present analysis. We briefly describe the minimal 
model, we discuss the front velocity and the coupling be- 
tween the front and the non-ionized exterior region, and 
we cite several properties of planar uniformly translating 
fronts required for our further analysis. 



A. The minimal streamer model 

The model for negative streamers in simple non- 
attaching gases like nitrogen and argon as used, e.g., 
in 0, [^, Ho, [O, E, Ej [iJ, ^M consists of a set of three cou- 
pled partial differential equations for the electron density 
a, the ion density p and the electric field E. In dimen- 
sionless units the model reads 



5tcr- V-(crE) = cr|E|a(|E|), 

dtp = cr|E| a(|E|), 

V-E = p-o-, E = 



-V(t>. 



(9) 
(10) 
(11) 



A general discussion about dimensions for this model can 
be found, e.g., in 0, SHll^- The first two equations are 
the continuity equations for the electrons and the ions 
while the last is the Coulomb equation for the electric 
field generated by the space charge p — aoi electrons and 
ions. Here the electrostatic approximation for the electric 
potential was introduced [y]. Notice that the electron 
particle current is here taken as the drift current only 
— (tE, neglecting electron diffusion. For the effect of a 
diffusive contribution, see the discussion in the introduc- 
tion as well as [43|. The ion particle current is neglected 
because due to their larger mass the ions are much less 
mobile. 

The term cr|E|a(|E|) is the generation rate of addi- 
tional electron ion-pairs; it is the product of the absolute 
value of the drift current times the effective cross section 
a taken as field dependent; a commonly used form is the 
Townsend approximation 



a(lE| 



^e-i/|E| 



(12) 



However, our analysis holds for the more general case of 

«(0) = 0, ^>0, 
ax 

(13) 



a(x) = l-O' ^ 

a;»l \ X 



B. Ionization fronts and their velocity 

We consider planar ionization fronts propagating in 
the positive x-direction into a medium that is completely 



non-ionized beyond a certain point Xf{t) 

a = Q ~ p for X > Xf{t), 

a > for X < Xf{t), (14) 

and we work in the comoving coordinate 

C = x~Xf{t). (15) 

As is easily verified by rewriting Eq. ([9]) in the comoving 
coordinate ^ with the help of pT|) as 



dt 



&- 



E • Vcr = cr (p-(T+|E| a(|E|)) 



(16) 

where x is the unit vector in the x-direction, the velocity 
of the front is determined purely by the value of the field 
E^ on the discontinuity of the electron density 



Vf{t)^^ = -E{xf{t),t)^E\ 



(17) 



Notice that, in order to simplify notations later, E~^ is a 
positive quantity. Obviously, this relation between front 
velocity and local field is not restricted to planar fronts 
only, but holds generally: 



v/(t) 



dt 



-E(r/(i),t) 



(18) 



As the space charge p — a never diverges, the field E is 
continuous across r/. 

C. Coupling the front to the outer non- ionized 
region 

In the purely one-dimensional planar setting, the field 
does not vary ahead of the front: dxE = ioi x > Xf{t). 
On the other hand, the curved ionization front around 
the tip of a streamer finger creates the characteristic field 
enhancement ahead of the streamer tip. However, the 
region x > Xf{t) can be left purely to the outer scale 
analysis of the Laplace equation 



V 



0, E = -V0 for ^ > 



(19) 



and has no further influence on the front solution ex- 
cept for determining the field E+ at Vf (t) . The matching 
between inner and outer region here is completely con- 
centrated on the line r/(i). As also remarked already in 
the introduction, these two spatial regions only decouple 
if the electron diffusion D (that is included in numerical 
streamer models) is neglected. Putting D — eliminates 
the leading edge and the pulled dynamics of the front 

am. 



D. Uniformly translating planar ionization fronts 

If _E+ is time independent, dtE~^ = 0, the front propa- 
gates uniformly with velocity v = E^ . Uniformly trans- 
lating front solutions of the minimal model (l5|l- (fTT1) in 



a field E^ depend only on the comoving coordinate 
S^ = X — vt. The y h ave been discussed in many previ- 
ous papers p, [y, H^], for a recent thorough discussion, 
we refer, in particular, to section II of [33] and for the 
comparison of fronts with or without electron diffusion 
D to section 2 of [Hi. 

Analytically the front solution in the ionized region 
^ < is given implicitly by the equations 



a{E) 



E^ 



E+ \E\ 



P{E), 



p{E) = / dpa{n), 

J\E\ 

dpE+ 



< 



-E+ 



ii) Ai p(Ai) 



(20) 
(21) 
(22) 



Fig. [^ shows the shape of the solution for E~^ — 1. At 
the shock front, the electron density jumps from zero to 



lim cr(£) ^ E+aiE+), 



(23) 



while ion density and electric field are continuous across 
the interface 

piO^OiO, E{0 = -E+ + O{0 for e<0. (24) 
For ^ ^ — oo the densities approach the value 



= cr(— cx)) = / dpa{ii), 
Jo 

= p{-Oo) = cr". 



(25) 
(26) 



which means that the space charge density p—cr vanishes, 
and the electric field vanishes as well 



E- = E{-oo) = 0. 



(27) 



Far behind the front, where E is so small that a(|E|) 
0, the profiles of E and a decay exponentially: 



e^^ + 0{e^^'^), (28) 




where 



V 



1 

E+ 



— - / dna{p) 



E+- 



(29) 



For more details on the asymptotic behavior at ^ = 
and ^ ^ — oo, we refer to section II of 42]. 

The exponential decay in ([28|) suggests a natural defi- 
nition of the front width as 



V 



(30) 



The length scales within the front are further elaborated 
in section 6.3.2 of |43| where next to l~ also another scale 
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FIG. 2: Densities and field for a uniformly translating planar 
ionization front in a far field \Eoo\ = E'^ = 1 as a function of 
the comoving coordinate ^ (I15() . The front moves to the right 
with velocity v = E^ . Upper panel: electron density a (solid 
line) and ion density p (dashed line); lower panel: electric field 
(solid line, axis on the left) and space charge density p — a 
(dashed-dot line, axis on the right). For a(|Ej) the Townsend 
approximation (|12|l was used. 



E. The minimal regularized boundary model 

The movin g b oundary model derived and evaluated 
in [H m [11,1331 is given by Eqs. ©-(El) together with 
the ideal conductivity approximation (/> = const, in f2~ 
that will further be discussed in Sec. |Vl This model 
can be considered as a minimal regularized boundary 
model. While the model with the boundary condition 
0+ — 0^ = treated in [30| lacks regularization and for 
generic initial conditions breaks down within infinitesi- 
mal time, the model of [32, j33, ^, JiTj according to our 
analysis does stay regularized as long as the initial con- 
tour is infinitely many times differentiable. It consists of 
the exact description of the non- ionized region p9| . of 
the exact relation between velocity and local field p^ 
and furthermore of two approximations for the full dy- 
namics of the PDE's (in])-((ni)- First, the jump of the 
electric potential across the front is modeled using the ap- 
proximation p3p of a planar uniformly translating front 
in the limit of large field E+; corrections to this behav- 
ior will be studied through the solvability analysis from 
section IIIII on. Second, the streamer interior is modeled 
as ideally conducting: — const, in fi^. This approxi- 
mation will be discussed and improved in Sec. |Vl 



t^ — l/a(_E+) was defined that characterizes the inner 
front region close to ^ = for D = 0. (Nonvanishing 
electron diffusion D introduces another scale ir,^ The 
length scale t^ is obviously of no relevance for matching 
of the outer scale, but t^ is the appropriate quantity. 

From the uniformly translating planar front solution, 
the boundary condition ([3]) was derived in [33] as follows. 
Integration of the electric potential, d(^(j){C) — ~-E'(C)> 
from — (X) to leads to 



0(0) -0(- 



dCEiC). (31) 



Substituting the variable -E(C) = —x with the help 
of (im, and using the limit (|27p for the lower limit of 
integration, we obtain 



(b+~r = Qo{E+), Qo{E+) 



dx ■ 



E^ 



p{x) 



-, (32) 



where the function p{x) is defined in ([2T|) . As discussed 
in detail in [33|, the function Qo depends nonlincarly on 
-E+ (see Fig. [5]), but asymptotically 



E+ for E+ > 1. 



(33) 



Finally, let us also recall an important property of the 
planar front solution which is a useful relation used later 
in the paper 



di 



aiOE{O=0. 



(34) 



This relation is obtained by taking the derivative of 
Eq. (122D and using Eq. ^. 



III. BOUNDARY APPROXIMATION FOR 
SLIGHTLY CURVED FRONTS 

To derive curvature corrections to the minimal regular- 
ized boundary model (section III Ep from the underlying 
PDE's described in section lTl Al we follow the general ap- 
proach as described, e.g., in 21. 2^[2^.[27l.[28l.T41. 45] for 
solidification fronts in supercooled melts, nonlinear dif- 
fusion fronts or interface between two inmiscible fluids, 
and adapt it to our problem as necessary. 



A. Perturbative expansion about planar fronts 

In a first step, we expand slightly curved fronts about 
the uniformly translating planar fronts that were recalled 
in section Hi Dl Weak curvature means that the width of 
the front £~ (|30p is much smaller than the mean radius 
of curvature R = 1/k of the front. The quantity 



k£: 



(35) 



is then a small parameter. The inner region of the curved 
front can be expanded about the planar front by intro- 
ducing a curvilinear orthogonal coordinate system, with 
coordinates (77, (^) locally tangential and normal to the 
moving front. We expand the inner region as 



v{f],t) 



(t)Q{C) + e(j)i 
cro(C) + ecri - 
Po(C) +epi - 
vq + evi + ■ ■ 



(36) 



where [(/'o(C)> fo(C)i Po(C)] and vq = E^ are the planar 
uniformly translating front solutions (HO])- (US) within a 
fixed field E+ . 

Assuming that the variations in the transversal direc- 
tion 77 are of the order of the radius of curvature i?, it is 
useful to rescale the transversal coordinate as 77 = erj. A 
divergence V • J then can be expanded in e as 



(37) 







V- J = 


9 T 9 T 


and the 


Laplac 


ian V'^cf) 


as well: 




V^ = 


92 


d 2 d^ 




= 




4 + 0(.2), 



A 



4 



(38) 



Therefore, the equations for [0i,(Ti,pi] and vi depend 
only on C as common. 

Substituting the expansion ()36|) in the minimal model 
equations ©-(HI]) we find to 0{e'^) the planar front 
Eqs. dlOll-llll). We obtain to 0(e) the equations: 

(po - 2o-o + /o) cTi +(Tq pi + (di^aQ - (To/o)£^i 

+ {vo + Eq) Oqcji = -vidi^ao, (39) 

/o cTi - aofoEi + vq dc^pi = -vid^pa, (40) 

-(Ji+pi-dQEi=\E^, (41) 

El = -ac0i, (42) 

where we used 



J{\E^ + Ei\) = U-Eifi 



(43) 



and where 
/o 



J{\Eo\) = \Eo\aa = \Ea\a{\Eo\), (44) 

dJ{x)U=\E,\ = ao(l - l/E^). (45) 

These equations are linear for {ai, pi, Ei,(j)i) with inho- 
mogeneous terms vidt^ag, vidt^po and Eq. The structure 
of this system of ODEs reads 

C-Vi=JV, £ = M--Dd(, (46) 

j^_ I -vid^po 



Ui 




AE'o 


/ Po - 2(70 + /o ctq a^^CTo - cTo/o \ 



(47) 



M(C) = 



/o 



-'^/^ 







D(C) 



V 0-1 

/-(7;o + £^o) 0\ 

-Wo 
10 



(48) 



0/ 







(49) 



1/ 



B. The structure of the solvabiHty analysis for 
streamer fronts 



In this subsection, using the zero modes of the adjoint 
operator £* of £, we formally derive corrections to the 
boundary relations derived from the planar front solu- 
tion. The explicit expressions of the coefficients in these 
relations are calculated in later subsections. 



1. Formal procedure of solvability analysis 

The boundary conditions are evaluated through a solv- 
ability analysis. We first introduce its general formalism. 
Consider a system of inhomogeneous linear ODEs 



£ • Ui = A/", 



(50) 



where £ is a linear differential operator that can be writ- 
ten as 



£ = M - D^c, 



(51) 



and where M and D are matrices whose entries depend 
on C,. The adjoint operator of C is given by 



C* =M* + a<;(D*) + D*9(;, 



(52) 



where * stands for the transpose of the matrix. Assume 
that a zero mode, U* , of C* is known 



£* ■ U* = on a<C<b. 



(53) 



Multiplying Eq. ([50]) with (U*)*, integrating over [a,b] 
and integrating by part, one easily verifies that 



dC (U*)*A/'=- (U*)*DUi 



(54) 



Here the arbitrary integration boundaries a and b were 
introduced to control possible divergences of the integral 
and the boundary term. The relation ([54)1 is the basis of 
the solvability analysis performed in this work. As we will 
show in the next sections, taking the limits a — > —00 and 
6 — > 0, relation (|54p yields linear relations between the 
perturbed electrostatic field and potential on both ends 
of the inner region, thus providing boundary conditions 
for the Laplace equations @ and (O in the outer regions. 
It should be remarked that the different coupling to the 
outer regions on both sides of the front is reflected by the 
asymmetric integration boundaries. At C = 0, there is an 
abrupt transition from the inner front region to the outer 
non-ionized region as described in section Hi CI here the 
matching region shrinks to a line. For C — > —00, there is 
an exponential decay to charge neutrality and the outer 
ionized region as discussed later in section FVl 



2. Zero modes of C and MBC equations 



The asymptotic behavior of zero modes at C, ~^ — oo 



Denoting the components of the zero modes as 



U* 




(55) 



relation ((54|) . with a = Cc and & = 0, leads to the relation 

(56) 



[(wo + EQ)ai^a- + vopiipp - iPeEi 
~Avi + BX, 



1° 



where 

A ^ I dC[i;adccro+i^pdcpo], (57) 

"'Co 

B = f dCEo^'E- (58) 

"'Cc 

In order to obtain the MBC equations we have to con- 
sider the limit Cc ~* — oo. This limit has to be analyzed 
with caution due to the mutual cancellation of diverg- 
ing terms. In order to extract from Eq. (j56p appropriate 
MBC equations, we use some properties of the perturbed 
field and densities. First, from Eq. (fT8|) and the expan- 
sion (1551) we have 



vi 



--^i(0). 



(59) 



Second, asymptotic analysis of Eq. ([55]) for ( -^ — (X) 
leads to 



(Tl(-Oo) = pi(-Oo), 



(60) 



which ensures charge neutrality at C = — oo to order e. 
Third, we show below that ijJaiO) does not diverge. Since 
from Eq. ([18]), vq + Eq vanishes at C = 0, the finite value 
of i7i(0) does not enter the MBC equations. Finally, since 
p is a continuous function and po vanishes for C > (|26p , 
Pi{0) must vanish as well: 



Pi(0)=0. 



(61) 



With these relations, Eq. 



takes a simpler form 



- [A + Vb(0)]«i +BX = EiiQ^EiCc) + MCc)MCc) 

- ^4O)0i(O)-«oai(Cc)[V'.(Cc)+V^p(Cc)], (62) 

where we assumed |Cc| 3> ^q, such that £'o(Cc), cro(Cc) 
and po(Cc) reach their plateau values (pS)) . 

In order to extract the actual MBC from Eq. ((62)) . we 
must find the zero modes of £*. The number of inde- 
pendent MBC equations is related to the number of zero 
modes that give rise to converging integrals in the expres- 
sions for the coefficients A,B (see Eqs. ([57)) and ([58)) ). 
The calculation of these zero modes and the consequent 
evaluation of the required integrals are the subject of this 
section and the next one. 



The construction of the adjoint operator C* from 
Eq. ()52p is straightforward. The system of equations for 
the zero modes, Eq. ([55)) . is: 

(.fo - Cro)'lJja + Mp -'ipE - {vo + Eo)d^1pcr = 0,(63) 

aoipa- + -ipE - vod^ipp = 0,(64) 

{di^ao - CTo/o)V'<T - ctq/oV'p - V'^ + ^c^s = 0,(65) 

dc^^ = 0.(66) 

Equation ()66p immediately gives 

4'(p = Ccj, = const. (67) 

As we explained above, the zero modes which are relevant 
for the derivation of MBC equations are those that give 
rise to converging integrals in Eqs. ([57)) and ([55)) . It is 
thus essential to analyze first their asymptotic behavior 
in the limit C ^ — cxo. In this limit, (Jq = po = o^ and 
Eq = fo = fo ^ 0, and the system of Eqs. (|63)) - (|66)) 
reduces to 

- a'Tpa- -TpE - vodi^-tpo- = 0, (68) 

cr"V<T +ipE - vodc'^p = 0, (69) 

OqiPe = C0. (70) 

Equation ([70)) immediately gives: 

'ipE = C^C + CE , (71) 

where ce is another constant. Substituting Eq. (|7ip in 
and integrating we obtain 



i>a 



CE 



c<f,vo 
(a-)2 



3.C 



(72) 



where Ca is a third constant of integration. Finally, the 
asymptotic expression for ^p is obtained by substituting 
Eqs. ([71|) and ([72]) into ([69]) and integrating: 



_ CT /■ C(j) 

ipp = -c^e "0 + -^(^ 



-p, 



(73) 



where the fourth constant of integration Cp was in- 
troduced. The set of four independent constants 
C(j, Cp, cb, C0 corresponds to four independent solutions 
of the set of four ODE's ([n3])-([M])- 

The exponential divergence of the coefficient of c^ in 
equation (|72p and the exponentially decaying profile of 
(To(C) with the same length scale (|28)) imply that the 
integral ([57]) converges only if Ca = 0. This condition 
reduces the number of relevant zero modes, and hence 
the number of MBC equations derived from Eq. ([S^ , to 
three. 



4- Initial conditions for zero modes at ( -^ — oo 

In order to extract the three relevant zero modes 
from Eqs. ([S5)) - ([5^ in the whole interval ( € 
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(— cxD,0), it is natural to specify initial conditions 

[lpa{Cini),tpp{Cini),'lpE{Cini),'4'4,{Cini)] at SOmC ( = (ini < 

—£~ , such that the zero modes are properly described by 
their asymptotic behavior, Eqs. ^^, ((TH) and ^^. Ob- 
viously, our final result must not depend on our choice 
of Cini- Using Eqs. ^, ^ and ^, and imposing the 
condition c^ ~ we obtain the relation: 



Co-e 



-Cini 



■0(7 (Ch 



V'g(Cini) _ 1p4,iCmi)vo 



From Eqs. 

V'.(C) 
V'p(C) 

0£(C) 



m 



0. 

(74) 



73|) and the constraint (f74|) we obtain 



^p(Ci„i) + ^^^^^(C-Ci 



-■0(7 (Ch 



Siniji 

i), 
^0(Cim)(C-Cir 



(75) 
(76) 

^)(77) 
cr 

From Eqs. ((75|) - ((77)) one can easily see that if 0,^(Cini) 
is nonvanishing, the zero modes diverge linearly with Q 
as C ~* — cx). Since at least one of the three required 
independent relevant modes must have ^'^(Cini) 7^ 0, we 
must include the computation of such a diverging mode 
in our analysis. For computational purposes, it is con- 
venient to choose the three independent orthogonal sets 
of initial conditions at C = Cini so that only one of them 
contains a nonvanishing value of ^-!,f,{C,ini)- It can easily 
be shown that choosing another set of three independent 
initial conditions, lead to identical results for the MBC. 
With our convenient choice we can immediately extract 
from Eqs. ((75|l - ((77|) the initial conditions at C = Cini (tak- 
ing Eq. ([7^ into account) for the two nondiverging modes 
[V'i*\V^^\^«,^«], where z = 1,2: 



and 



/V'^'^(Cini)\ 
V'^'^(Cini) 
V'i'^(Cini) 

V#(Cini); 



/4'^(Cini)\ 

V'f(Cini) 
/,(2)/ 



(78) 



1 






(79) 



V'r(Cini) 
VV'i'^(Cini); 

where the linearity of Eqs. (p5)l - (|66p was used to scale 
one of the coefficients to unity in each initial condition 
(|78| and ((79| . The last set of initial conditions, for the 



linearly diverging zero mode [^o- tWp , we 1 0^ ] i i^ ^^" 
termined by requiring it to be orthogonal to the two vec- 
tors dZH]) and dlHl) (and to satisfy the constraint (fM)) ). 
This condition gives, up to a normalization factor: 



/4'^(Cini)\ / 



0f(C 

V-f^lCin 
,,(3) 



V4''(Cini)/ 



\ 



Vf7(i + (--ny 



Thus, in order to obtain the three zero modes relevant 
for the MBC equations, from Eq. (|62p . we have to solve 
Eqs. ((63|) - ([66|l in the regime C G (—00, 0) under the ini- 
tial conditions, Eqs. ([75]) - ([50)1 . Computation of the lin- 
early diverging mode [V'o- \i>p ■,'Ve ^"Vl ] and extract- 
ing nondiverging MBC equations from Eq. ([^^ corre- 
sponding to this mode involves nontrivial technical com- 
plications, which for the sake of clarity we defer to the 
following section. We first discuss the two MBC equa- 
tions corresponding to the two converging zero modes 



5. The nondiverging zero modes 

In this subsection we derive a general form of two MBC 
equations. The exact expressions for the coefficients are 
deferred to section UlI CI As we show below, the general 
form of the two MBC equations we derive here is suffi- 
cient to deduce an important property: The perturbed 
electrostatic field Ei{—oo) does not vanish in the inner 
region behind the front. This observation emphasizes the 
necessity to derive the third MBC equation with the aid 
of the linearly diverging zero mode. 

For the two nondiverging modes we can take Cini ~ 
Cc -> -00. Substituting Eqs. ^ and ^ in Eq. §^ 
we obtain the two MBC equations: 

^o'Ti (-00) = [^(1) + ^g^ (o))]„^ _ s(i) X , (81) 



and 



VqCTi 



(-00) + (T- ^1 (-00) = [^(2) ^ ^(2) (Q)]^^ - S(2) A . 

(82) 



Let us reconsider now our perturbative expansion (j36p . 
As we emphasized in previous sections, the only physi- 
cal source of perturbations to the planar front solution 
in our analysis are the transverse perturbations which 
correspond to nonvanishing k = Xe. Thus, we formally 
write the perturbation to front velocity as e wi , requiring 
it to vanish for e = 0. The calculations presented so far, 
however, do not involve any assumption on the actual 
value of vi, except evi <C vq. By keeping k = in our 
former analysis, we can thus consider the formal limit 
process: e ^ with evi small (compared to vq), but fi- 
nite. This simply means that we consider a perturbation 
of the electric field for a pure planar front which implies 
Ei{—oo) = 0. Repeating our former analysis with k = 
in Eq. (j38p we observe that the inhomogeneous term in 
Eq. ([41]) vanishes, leading to B = in Eq. ^E^j. Substi- 
tuting B = and £;i(-oo) = in Eqs. ^ and ([MD we 
obtain the nontrivial identity: 



VqCTi 



(-00) = [^(1) +7AW(0)]«i = [A^'^ + ^p'^^\0)]v. 



(80) which implies 



^W-fV'i'^(0)=A(2)_^_^(2)(Q) 



(83) 



(84) 



This important identity will be confirmed in Sec. IIIIDI 
after finding the exact solution for the converging zero 

modes, and evaluating the actual values of A^*^ +^^ (0), 
for i = 1, 2. 

Furthermore, substituting the identity (|84p in Eqs. 
([8T|l and (l82|) we obtain the equation: 



5(2)_5(1) 

Ei{-oo) = A 



(85) 



This relation implies that the electric field behind our 
curved fronts does not vanish if i?^^^ ^ B^^\ 

Finally, let us use the relation ((85|) to derive a cor- 
rection to the boundary condition ([3]). Integrating the 
equation 9(;0i = —Ei, we get: 



0i(Cc)-'^i(o) 



dCEiiC). 



(86) 



Since Eq. (|85|) implies that £'i(C) reaches a constant 
value Ei(—oo) as C — » — oo, substituting this in Eq. 
we obtain the asymptotic form of (f>i{C)'- 






(87) 



where Wi is a constant. This relation will be used in 
Sec HVn 



C. Exact solutions for the zero modes 

In this section we find exact solutions of Eqs. 
(|66| under the three sets of initial conditions, Eqs. (jTS]) - 
(IMl) , thus obtaining the three zero modes relevant for the 
MBC equations derived from Eq. (p^ . The calculations 
are lengthy and many technical details are delegated to 
appendix [X] Since the final expressions we obtain are 
considerably complicated, we support our analytic cal- 
culations by comparison to numerical solutions of the 
system of ODEs (|63|) -(|66 l) . To avoid cumbersome no- 
tations, we omit the superscript distinguishing between 
different zero modes in general expressions that apply for 
all modes. 



1. Analytic solution 

Our starting point is Eqs. (p5)l - (|66p . With some al- 
gebraic manipulations (see appendix E]) . these equations 
can be brought to the form: 



^vo + Eo)d^^{^^ + ^p) + fodcii^^ + i,p) 

C^Eq 



vo 



= 0, 



(88) 



Eodc^a = -vod^i^, + Vp) + h{^. + ^p) , (89) 

^cV'E = 0-o/o(V'a + i>p) - {dcao)i)a + C0 . (90) 



Notice the structure of this set of equations: Eq. 
is an inhomogeneous first order ODE for ^^(V'o- + V'p) 
and can thus be solved analytically. Once this solution is 
substituted in Eq. ([55]) . this last equation becomes also 
an inhomogeneous first order ODE for ip^, and is thus 
also analytically solvable. Finally, substitution of the so- 
lutions to Eqs. ([55)1 and (|5^ yields an inhomogeneous 
first order ODE for ipE, which is again analytically solv- 
able. In the rest of this section we follow this procedure 
to obtain general analytic expressions with appropriate 
constants of integrations, which are then determined by 
imposing the initial conditions, Eqs. ([75)1 - ([50]) . 
The general solution of Eq. ([55]) has the form: 



5c(V'.+^p)(C) = C(C)e«('^); ff(C) 



C 



dx- 



foix) 



(91) 
where the expression of C(C) is given by (jA12[) . Inte- 
grating from (^ = — oo to C (the sum ip^ + ipp does not 
diverge at C = — oo, see Eqs. ([75]) and ([76|l ) and using the 
change of variables y = Eq{x) to calculate the resulting 
integrals, we obtain the final form 



(V'a + ^p)(C) 



^ap 



vo Jo p{xy 



where 



Cap = (V'ct -l-'0p)(-oo). 



(93) 



The explicit steps leading to this expression are given in 
appendix lA II 

Substituting Eq. ^^ in Eq. ^^, we obtain after 
integration and using the same variable transformation 
y — Eq{x), the complicated expression (|A29[) given in 
appendix lA 21 This expression can be written in the fol- 
lowing compact form 



MO = ^Cini 



aaiO+baiOe 



ACi, 



(94) 



where a^iC) is given by Eq. (jA32[) and ^^(C) is a regular 
function of C whose exact form does not affect our analy- 
sis since we will consider later the limit |<^i„i| ^ 1 and the 
last term of Eq. ([M)) will be as small as we want. Notice 
that if C0 7^ 0, as is the case for the linearly diverging 

mode [tpa , V'p , V's 5 ''/'i ] ) the first term on the RHS of 
Eq. ([M)) diverges as the auxiliary parameter C^ini — > — oo. 
As we show below, this singular dependence on Cini can- 
cels out and does not appear in the final form of the MBC 
equation extracted from Eq. (|62[) corresponding to this 
mode. Since in deriving Eq. ([5^ in Sec. IIIIB 21 we as- 
sumed, however, that ipa-{0) is finite, we must carry the 
(^ini-dependent terms that appear in Eq. (|94p and other 
related expressions throughout our analysis. 

Substituting Eqs. ^ and ^ (whose fuU form is 
given by Eq. (|A29p ) in Eq. ([501 . integrating and using 
again the variable transformation y = Eq{x) we obtain 
the cumbersome expression (IA4ip given in appendix IA3I 
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(notice that to simplify this expression, we used for the 
first time the exphcit form (fT^ of a{x)). This expres- 
sion can also be written in the following compact form 
similarly to Eq. ([M|) 



^EiC) 



a 



+ aB(C)+&£(C)e^^'"' , (95) 



where asiC^ is given by Eq. (IA44P and 6b (C) is a reg- 
ular function of ^ whose exact form does not affect our 
analysis. 



2,. Comparison with numerical results 



To check the analytic formulas ((9^ . ((M|) and 
with the complicated expressions, Eqs. (jA24[) . (|A29p and 
(|MT|) . we used MATHEMATICA to obtain numerical 
solutions of the zero modes Eqs. ((63)) - ([66l) under the ini- 
tial conditions (l78 |) -([80 |) . for few representative values of 
Cini- In Fig. [31 we present the numerical solution obtained 
with the initial conditions for one of the two nondiverging 
modes, Eq. (|79)) . evaluated at Cini = —40 for Eao ~ —2. 
In Fig. m we report the numerical solution for the lin- 



1.5 



1.0 



0.5 



-1.0 



-1.5 



-2.0 




-16 -14 -12 -10 



-8 -6 



FIG. 3: Comparison between numerical solution of Eqs. (|63p - 
(|66|) and the exact solution reported in Sec. lIIICl for the initial 
condition ((791), with Eoo = -2 and Cini = -40. 

early diverging mode, obtained with the initial condi- 
tions, Eq. ([501) , evaluated at Q^i = — 50 for Eaa = — 1- 
Figs. [21 and m exhibit a excellent agreement between our 
analytic expressions and the numerical solutions. 




FIG. 4: Comparison between numerical solution of Eqs. (|63|l - 
(|66p and the exact solution reported in Sec. lIIICl for the initial 
condition ([80l) . with Eao = — 1 and Cini = —50. 



these expressions are given in appendix \^ Using the 
definition of A, Eq. ([57|l . and the exact expressions of 
-0CT + V'p, Eq. ([5^ . and of V'a, Eq. ^ (given in its full 
form in Eq. (|A29p ). we obtain the explicit form for A, 
Eq. (IA57P (see apeendix IA4p . This expression can be 
written in the compact form 



_ (7O(O)C0 
-^ — ~ S,ii 



+ aA + bAe^'^"^\ 



(96) 



where ua is given by Eq. (|A59p and the exact form of 
hA is not required for our analysis. The exact form of 
the coefficient A + 4'e{^) is given in Eq. (|A60p which 
is obtained using the expressions ([Ml and ([551 for A 
and i}E respectively (see Eqs. (|A4ip and (|A57p for their 
complete expressions). The compact form for A + V'b(O) 
is given by 



A + Vb(0) = aA + aE{0) + bE{0)e^'^" 



(97) 



Notice that this last expression does not diverge when 
Cini -^ — oo. The term qa + asiO) is given by Eq. (|A62p . 
Using the definition of B, Eq. ([551) . ^nd the expression 
([Ml for Ve(C)j we obtain the comphcated form of B, 
Eq. (|A68p . The compact form for B is given by 



S = ^Cini + aB+6Be^'^" 



(98) 



where qb is given by Eq. (|A70p and the exact form of bs 
is again not required our analysis. 



3. General expression for the relevant coefficients 

In this section, we use the analytic solutions derived 
in Sec lIIICTl to provide compact expressions for the co- 
efficients A, A + V'b(O) and B, required for deriving the 
MBC equations from Eq. ([0^ . The complete forms of 



D. Completing inner region analysis 

We here complete the inner region analysis by return- 
ing to Eq. ((nH and using the results obtained in Sec. lIIICl 
to derive MBC equations for the curved front. 
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First, let us rewrite Eq. ((62)) in the form 

- {A + ^EmVl +B\ = EiiCc) i-MCiniW' 

+ c^(Cc - Cini + ^)) + c^MCc) - MO)) 

a J 

- Co-pWoCri(Cc), 



(99) 



where we used Eqs. (|67l) . (|77l) and ((93)) (taking aheady 
the Hmit Cc -^ — cxo). To further facihtate our analysis, let 
us first write the expressions of the coefficients A + -0£;(O) 
and B in the following form 

A + ?/>b(0) = 7lC0 +72Cctp +73V'a(Cini) 

+ 74^B(Cini) + 6B(0)e^';'-, (100) 

B = i-^Cini+(3l)c^+fhc„p+p3^PaiCini) 



+ p4i'E{Cini) + bse^'^" 



(101) 



The coefficients "fj and Pj, j = I,--- ,4, depend only 
on _E+ and can be easily extracted from Eqs. ((A62p and 
(EZOj: 



7i 



E+ 



I (wq - a;) 



(102) 



dx[a(x) -«(£;+)]- 
p(a;)^ 

72 = i'oa(£;+) (103) 

73 = o~ (104) 

74 = 1. (105) 



/3i 



I'D 



E+ 



dx ("°-"^"(") rdy^^^ (106) 



rfa; 



a:p(x) 



/32 



cr ./o a;p(a;)2 

(wo - a:)a(a:) 



P(2/)^ 



[P(2^)W0 - CT (Uo - 2^)] 



B+ 







dx- 



E+ 



xp{x) 



vo- X 



Ps. = v^-cj I dx- 

Jq P[x) 

.E+ 

f3i = - dx " ■ 

Jo P(x) 



(107) 
(108) 
(109) 



The contributions of last the term on the RHS of 
Eqs. ()100p and ()10ip can be made arbitrarily small by 
choosing K'inil sufficiently large and can thus be neglected 
in the rest of our analysis. 

The first MBC equation is obtained by using the val- 
ues of the nondiverging zero mode corresponding to the 
initial condition ([75]) in Eqs. ((Mj), pTO)) and ((TUT)) . We 
obtain 



Wocri(-oo) = 72^1 - /32A, 



(110) 



where we took the limit (c -^ — cxo. The second MBC 
equation is obtained in a similar way by using the values 
of the second nondiverging zero mode corresponding to 



the initial condition ([79]) in Eqs. ([99|), ((TOO)) and pOT)) . 
Using Eq. (jllOp . we obtain 



a 



(111) 



As was pointed out already in Eq. ((55)) . this last equation 
implies that the electrostatic field does not vanish behind 
the front. Notice also that, as mentioned in the discussion 
following Eq. ((84)) . it is now easy to verify explicitly the 
identity 



^(1) + V^i,') (0)^^(2) + 



i>%\o) 



12- 



(112) 



Using the nondiverging zero modes to extract MBC 
equations ()110|) and ()111|) from ((99)) was relatively 
straightforward since all terms are finite at both end 
points C = and ^ = Cc — > — oo- Naively, one may 
anticipate that since ^£; (Cc) ^ Cc it is impossible to ex- 
tract physically meaningful MBC equation correspond- 
ing to this mode from ()99p . We show however below that 
taking appropriately the asymptotic limits C,c -^ —00 and 
the auxiliary parameter \C,irii\ S> 1 leads to cancellation of 
the diverging terms and to an additional nontrivial MBC 
equation. 

In order to derive this third MBC equation, using the 
initial condition ((50)) . we first substitute ((57)) in Eq. 
and use Eqs. ((T00| and ((TUT)) to obtain 



[7lC0 + cr (72 +73) +74] Wl 
„2 
( — Cini + /3l)c0 -t- CJ-{[i2 + /33) + (ii 



= c^Wi-a vaai{Cc) + Ei{Cc)[c4,vo/a ~ {a )^] 

- C^Cini^l(Cc), (113) 

where c^ is given by ((50)) (using Eq. ((07)) ). Notice that 
the divergence of ip^ and ipi with Cc cancels out provided 
Cc| is large enough such as £'i(Cc) has reached its plateau 
value. We can now take the limit Cc ~^ ~oo such as 
the limits in Eqs. (jllOp and piip are reached and can 
be used in Eq. (|113p . We then also notice that the two 
terms proportional to Cini in Eq. ((113p cancel. After some 
manipulations, our third MBC equation can be written 
as 



Wi 



-UJ\ V\ + W2A 



where 



wi = 71 + 



"0 



r. Vo f vo-x v^ 

^2 = Pi dx + _ 

o- Jo pyx) [a Y 



(114) 

(115) 
(116) 



In Fig.[5l we present the dependence of the coefficients 
in Eqs. ^, ITTT)) and ((TTi)) on the variable E+ . The 
relations between these functions and the functions Qi, 
used in Eqs. p26p and (|127p . are derived later in this 
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FIG. 5: vllo- = Q2{E+), a;i, lu2 = -Qi{E+) and Qo(£+) 
as a function of vo = E^ . See Sec. U for the relation between 
the coefficients Qi and the MBC equations. 
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FIG. 6: The coefficients 72(-B^) and /?2(£+) which determine 
the modified electron and ion density behind the front (see 
Eq. dnOj) 



section. In Fig. [HI we present the dependence of the co- 
efficients in Eq. (|110p on the variable E'^ . 

Notice the structure of the MBC equations pTU|) . (fTTTI) 
and (|114p derived above. As we will show below, the last 
two equations give rise to the MBC equations (|126p and 
(|127p . whereas the first one determines the excess density 
of ionized matter behind the curved front which is thus 
"slaved" by the dynamics specified by Eqs. (fT24|) - (fT28)l . 
(Hm, dnU and (UHl), in view 
143 1 for planar front solution. 



Let us discuss Eqs. 
of results obtained in Ref. 
For purely planar front, the only possible perturbation 
of the quantities (ctq: Po and E^) are due to modification 
of the electric field ahead of the front -Eo(O) = E^o- As 



already discussed in Sec. IIIIBSi in our formalism this 
is equivalent to keeping in Eqs. (jllOp and (|114p only 
the terms proportional to vi. In this case we also have 
-Ei(— oo) = 0, which implies that Eq. ([57)1 reduces to 

0i(-«))-0i(O) = I^i. (117) 

A first test concerns our Eq. (jllOp . In Ref. [43| the 
following relation was derived 



da 
dE+ 



dao{~oo) 
dEoiO) 



-a{E+), 



(118) 



with -Eo(O) = —E^. This relation actually follows di- 
rectly from the definition of a^ (|25| . Since for infinites- 
imal perturbations d(j{}[—oo) = (Ti(— oo) and dEoiO) — 
Ei{0) we obtain from Eq. (fTTS)) 



cri(-oo) = -a{E+)Ei{0) = a{E+)vi. 



(119) 



Keeping on the RHS of Eq. (jllOp only terms proportional 
to vi and using (|103p we recover Eq. (|119p . 

The second test stems from comparing our MBC 
Eq. (|114p to the planar front relation Eq. ([321) , which for 
infinitesimal perturbations, 4>i and Ei, recast the form 



M-^) - MO) = ' 

A simple computation leads to 
dQQ{E+) _ 1 f^^ dx 



dQo{E+ ) 
dE+ 



■Vi. 



(120) 



dE+ 



a{E^ 



p{x) 



<E+) 



E+ 



dx 



E^ 



p{xf 
(121) 



Keeping again only terms proportional to vi on the RHS 
of Eq. (|114p . using Eq. ()117p and comparing with (I120p 
implies that the following equality 



dQojE^ 
dE+ 



= ^1, 



(122) 



must be satisfied. To check this we use the definition 
(fTT5)) together with pI7^ to obtain 



Wi 



-a(uo 



dx 



dx 



Vo 



P{x)^ 

a{x){vo — x) 

p{xY 



a^ 



(123) 



Using Eq. ([2T|) . integrating by parts the second integral 
(using dx{l/ p{x)) = a(x)/p(xy') and taking the appro- 
priate limit of boundary terms, Eq. (jl23p becomes identi- 
cal, as required, to Eq. (|12ip . These last two tests demon- 
strate consistency of our calculation with the previous 
results of [H. 



IV. A PHYSICAL DISCUSSION OF THE 
CURVATURE CORRECTIONS 

A. The new moving boundary conditions 

In this section, we derive a refined version of the 
boundary model which includes curvature correction. We 
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first write the equations of the model and we discuss it 
below. It consists of Laplace equations for both outer 
regions 



V^f/) = in S1+ 
V2(/) = in SI" 

and the moving boundary conditions 



(124) 
(125) 



n.V0" = Q2{n-\/(j)+)K (126) 

(j)+-(t>- = go(n-V0+) + gi(n-V0+)K (127) 

vn = n-V0+, (128) 



where the coefficients Qi depend on the electrostatic field 
ahead the front, and arc given by analytic formulas de- 
rived from the planar front solution as follows 



Qoiy) 
Qi{y) 
Q2{y) 



, y-z 

dz——, 
P{z) 

-^2(2/), 
y' 



<^ (y) 



(129) 
(130) 
(131) 



and where p, a~ and L02 are given by (|2T|) . (125]) and (I116P 
respectively. 

From inner analysis performed in the previous sections, 
and as already discussed in Sec. IIIC| we know that the 
electron and ion densities vanish precisely at the discon- 
tinuity line of the front part of the ionization front where 
the matching with the outer region 17+ has to be done. 
Since there is no charge in fJ"*" the only consistent equa- 
tion for the electric potential in this region is given by 
(|124p . On the other hand, we know also that the space 
charge density vanishes at the back part of the inner re- 
gion (see Eqs. ^ and JSO])). Consequently Eq. (fT25l) 
for the outer region fi" is also consistent with the inner 
analysis. Another possible equation for the electric po- 
tential in Qr is given by V^(/) = a — p with the necessary 
boundary condition cr — p = on F" in order to be con- 
sistent with the inner analysis. However this equation 
introduces new unknowns, the electron and ion densities 
in ri~, for which additional evolution equations would be 
needed. In Sec. |Vl we discuss the constant conductivity 
approximation introduced by Eq. (|125p . 

The boundary equation (jl26[) is obtained directly from 
Eq. piip . Indeed, as e — > 0, the inner region reduces to 
the curve F (the sharp interface) and 



Ei{-oo) = n- V^j" 
On the other hand we know that 



vq 



E+ = h- V(t)^ 



(132) 



(133) 



Consequently we have in the first order in e using the 
expansion ^ and Eqs. ([T^ . P^ . ([T^ and P^ 



n. S/(j)- = n- V( 



o- (wo) 



K. 



(134) 



Since n- V<^^ = and e = k/A, we just recover Eq. (|llip . 
The boundary equation (|127p is obtained using 
Eqs. dHZ]), (fTTTjl and (fTT4| which reads 



01 (0) - 01 (Cc) = wi«i - W2 + ^Xc A. (135) 

V cr (uo) / 

The quantity Cc lies in the inner region and will coincide 
with 0" when e ^ whereas 0-^0+. Using the ex- 
pansion dSH]) together with Eqs. ([321), dM]), (fTUll . (fT22| . 
P^ and P^ we obtain 

0+-0- = 0+-(/)-+e(0+-0-) 

= Oo(S+) - eVFi = Q^{E+) + e{u:^vi - c^sA) 
= Oo(n-V0+) + Oi(n-V0+)K, (136) 



which is just the boundary condition p27p . 

Eq. (I127P contains a curvature correction to the "ki- 
netic undercooling" boundary condition ([3]). This cor- 
rection might suppress non-regular solutions [38] of the 
minimal regularized boundary model. 

Eq. p26p has required us to take the finite conductivity 
of the ionized region into account (Eq. (|125p ): it expresses 
the field behind the ionization front in terms of the field 
ahead of it, of the curvature and of the conductivity in 
the ionized region. 



B. The field inside a propagating streamer 



Eq. p26p has an interesting direct physical conse- 
quence for streamers whose radius and field enhancement 
evolve slowly during propagation within a comoving co- 
ordinate frame. Indeed, simulations of single streamers 
show that the field inside the streamer is constant in very 
good approximation (see for example Fig. 4 in [33|). Such 
a constant interior field is a common ingredient of phe- 
nomenological streamer modeling. Here we can calculate 
this field: on the streamer axis we find 



E 



-Q2(E+)kz. 



(137) 



That the field inside the streamer might approach some 
fixed constant value, is a common interpretation of exper- 
imental results. Here we have derived an explicit relation 
for this field; more precisely, it depends on the conduc- 
tivity inside the streamer, on the front velocity and on 
the curvature. This explicit prediction will be discussed 
further and tested on simulations elsewhere 14611. 



C. Large field limit 

Since the width of the front £„ (wq), see Eqs. ([^ and 
([5(I)l . is a monotonic decreasing function of vq, we expect 
that the sharp interface will be a better approximation 
for large fields. Moreover, as seen in Fig. [51 the coeffi- 
cients Qi involved in the boundary conditions (|126p and 
p27p become linear in this regime which simplifies the 
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model. In this section, we compute the asymptotic be- 
haviors of the coefficients Qi . 

From the definition (I129P of Qq{vo) we have 



dQo ,. vq- X 
= hm 



1 



dvo x^vo p{x) ct{vo) 
Consequently, we simply have 

Qo{vo) ^ Wo for Wo -^ oo 



1 for Wo — > oo. (138) 



(139) 



Using the same method, the limit of the constant con- 
ductivity behind the front (|25p takes the form 



a {vq) ^ Vq for Vq —+ OO, 



since 



da 
dvo 



a{vo) '^ 1 for Vq -^ oo. 



(140) 



(141) 



We then deduce that the asymptotic behavior of Q2 (|13ip 
is given by 



Q2{vo) ^ Vq for Wo ^ 00. 



(142) 



The computation of the asymptotic behavior of Qi is 
a bit more lengthly. The expression of this coefhcient is 
given by five terms. Let us compute the limit for each 
one and use the notation qi with i = 1,2,- •■ ,5 with 
Qi = X]i=i 1i- Using the results obtained above we easily 
find that 



qi 



{a-y 



'Vq for Wo — > 00. 



(72 = ^Qo(wo) ~ Wo for Wo -^ 00. 

(7 



(143) 



(144) 



The third term, actually contained in the expression of 
/3i (I106p . is given by 



93 = T / dx 

o- Jo 



(wo + a;)(wo - x) _ vq 



xp(x) 



= T'?3(wo)- 

a 



(145) 

The factor in front of the integral tends to 1 for large wo. 
The asymptotic behavior of the integral 93 is obtained as 
above: 



^*=2 1im "°-" 



dvQ x^vo p[x) ol{vq) 

We then find 

53 ^ — 2wo for Wo -^ 00. 
The fourth term is given by 

(wo + a:)(wo -xY 



2 for Wo "> 00. 

(146) 



f-vo 

qi^ I dx- 
Jo 



xp{xy 



(147) 



(148) 



To obtain the asymptotic behavior we also use the deriva- 
tive: 

dl^ = 2 1im ("°'!)'=2 1im ("° " ^^ 



dvo x-tvo p{xy x~*vo p(x)a{x) 

2 

2 for Wo -^ 00. (149) 



a(wo)2 



We thus obtain 



(74 ~ 2wo for Wo ^ 00. 



(150) 



Which means that 9i + 92 + 93 + <Z4 "^ ^'((wo)*') for large 
Wo. Since we know from numerics that the asymptotic 
behavior of Qi is linear, this behavior is encoded in the 
last term: 

r° (wo - x)a{x) r (wo - y) , , 

With the help of the following changes of variables x = 
Wo 9 and y = vqp we get 



4 f\ (1 ~ q)a{voq) f" (1 - p) 
95 = -Wo / dq r / dp— r^. 

Jo qpivaq) Jo p{vopr 

From the definition of p{x) (PT|) we have 



(152) 



p{voz)=vo dpa{vQp), (153) 

and from the definition of a{x) (J12p we know that 



lim a{vQz) = 9{z), 

Vo — 'OO 



(154) 



where d{z) is the Heaviside's function. We then get 

f\ In(l-g) TT^ , ,^,,. 

95 ~ Wo / dq = — —Wo for wo -^ 00. (155) 

Jo q 6 

Finally the asymptotic behavior of Qi{vo) is given by 



)i(wo) ~ — —Wo for Wo ^ 00. 
6 



(156) 



For large electric fields (£■+ > 4), the boundary condi- 
tions (fT^ - P^ reduce 

n. V0~ = Kh-V(j)+ (157) 

(?;.+ - 0" = (1 - tt^k/G) n ■ V0+ (158) 

w„ = n-V0+. (159) 

D. Algorithm 

Let us describe now an elementary algorithm for nu- 
merical solution of the front dynamics with the MBC 
Eqs. (I124p - (ll28p . For simplicity let us consider a front 
characterized by a curve {xf{y, t),y), and far field bound- 
ary conditions: (i) (j) = at an electrode {x = 0, y) behind 
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the moving front, and (ii) ^ = —Eao a,t x ~ L ^ £^. 
For simplicity, let us further assume an initial state of the 
front {xf{y,0),y), where Xf{y,0) = xq +Sf{y), xq >C 
and \Sf"{y)\ <C 1 and where ^~ <C xq ^ L. For such an 
initial condition, the initial expression for Q2 in Eq. (I126[) 
is approximately Q2(|£'oo|) (up to correction 0{Sf"{y))). 
The front dynamics is then prescribed by the following 
numerical algorithm: 



1. solve Laplace Eq. (|125p in the region fi with 
boundary conditions (i) at (x = 0, y) and (|126|) 

at (a: -^ {xf{y,t))-,y) . 

2. Use the solution of step 1 to evaluate the potential 

(j)-{x^ {xf{y,t))-,y). 

3. Substitute the potential (j)~ evaluated at step 2 in 
Eq. (fT27)) . and solve the Laplace Eq. (fT24| in the 
region Q'^ subject to the boundary conditions (ii) 
and p?71) . 

4. Use the solution of step 3 to evaluate the gradient 
of the potential n- V(j)^{x -^ x f {y , t))~^ , y) . 

5. Advance the curve (xf{y,t),y) in a normal direc- 
tion n (towards il'^) at a rate n- V(/)+ found in step 
4. 

6. Modify the functions Qi according to Eqs. (|129p - 
([T3T|) in the MBC Eqs. p26)) and p27|) according 
to the value of n • V0^ found in step 4 and return 
to step 1. 



V. THE CONSTANT CONDUCTIVITY 
APPROXIMATION IN THE IONIZED INTERIOR 

A. The degree of ionization behind a front 



Eqs. P5l) . ([26)1 show that the degree of ionization 

a^ ~ P— Jn dii a{fj,) behind a uniformly translating 
front depends on the electric field strength E^ . In the 
appendix of 47], it is shown that precisely the same de- 
gree of ionization a^ = a^{E^) from Eq. (pS)) is reached 
behind a planar front even if it does not propagate uni- 
formly, if only the electric field E^ is the same at the 
moment when the electron discontinuity passes the point 
of observation. Behind a slightly curved front, the degree 
of ionization will be approximately the same. 

As a consequence, the conductivity further behind the 
front varies as a function of evolution history, as the im- 
mobile ions with their positive charge act as a memory. 
E.g., when a streamer emerges from an avalanche in a 
homogeneous background field Eback, the conductivity 
in the evolving channel will increase from cr^ (| Eback |) 
immediately after the avalanche phase to tT~(|Eonh(OI) 
where Eonh(i) is the enhanced field at the tip of the 
streamer at time t. This increasing electron and ion 
density along the axis of the streamer finger is typically 
seen in simulations, e.g., in Fig. 1 of [485 fo'" ^^ emerging 



streamer or in Fig. 3(a) of [Tjl for a fully developed long 
streamer. 

When an electric field is applied to a body with spa- 
tially varying conductivity, in general, space charges can 
form in the interior, since charge conservation dtq+^-j — 
with q ^ p — a and j = crE can be rewritten in the form 



[dt + 0-) (p - 0-) + E • Vo- = 0, 



(160) 



in agreement with ((9|l- (fTT|) . i.e., the space charge p — u 
will decay on the Maxwell time scale 1/ct only if there 
is no conductivity gradient SIcr in the direction parallel 
to the electric field E. However, in practice such space 
charge effects are typically small if the evolution is slow 
and consequently |Vcr| is small. 



B. The approximation of ideal conductivity 

Lozansky and Firsov in their Russian textbook and in 
their short English article [29'| suggested to neglect the 
weak field in the ionized interior completely together with 
possible weak space charge effects, and to approximate 
the streamer interior as ideally conducting, i.e., to assume 
infinite conductivity and 



const in il 



(161) 



Inspection of simulations show that this is a valid ap- 
proximation for the weak interior fields in comparison to 
the strong fields in the exterior. Therefore the approxi- 
mation was adapted in the model |32l . l33l . l36l . |37| , and its 
validity was further tested on simulations in section IV 

offal. 



C. The approximation of constant conductivity 

A shortcoming of the ideal conductivity approximation 
is that it does not allow for a potential drop along the 
length of the streamer which is an important observable 
in experiments. In fact, the electric field along the axis 
of a long streamer is frequently assumed to be constant; 
it is an important point of scientific discussion and an 
ingredient of coarse grained models for complete streamer 
trees. 

The approximation of constant conductivity in the in- 
terior is justified if the width of the streamer and the 
field enhancement at its head vary relatively slowly dur- 
ing propagation as it is the case in many simulations. We 
therefore here adapt the approximation of constant con- 
ductivity (j~ in the ionized streamer interior. According 
to Eq. (|160p . space charges p — a then decay on the time 
scale 1/(7^. For a front propagating with velocity v, this 
corresponds to a spatial decay on the length ^~ = v / a^ 
(1501) . generalizing the explicit asymptotic behavior pS)) 
derived for uniformly translating planar fronts. Behind 
this layer, there are no space charges and the Laplace 
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equation approximates the ionized interior as well: 

V^^^O, E = -V0 in n-. (162) 

From a mathematical point of view, a finite conductiv- 
ity in the interior might lead to an additional regulariza- 
tion mechanism. We have a model with two Laplacian 
fields on both side of the interface. This type of model 
is called Muskat, or Muskat-Leibenzon, problem in hy- 
drodynamics and describes the evolution of the interface 
between two immiscible fluids in a porous medium or 
Hele-Shaw cell under applied pressure gradients or fluid 
injection/extraction. The so-called Hele-Shaw problem 
is thus the o ne-p hase version of the Muskat problem, see 
for example [11 |M HH HI ■ 

VI. CONCLUSION AND OUTLOOK 

We have derived the curvature corrections to a mov- 
ing boundary approximation of streamer ionization fronts 
within the minimal model. We also have introduced 
a constant conductivity approximation in the ionized 
streamer interior. The matching procedure between in- 
ner and outer region followed the standard procedure at 
the back end of the front. However, ahead of the front, 
the matching region shrunk to a line. The solvability 
analysis was based on integrating from — oo up to this 
line. It involved three zero modes of the adjoint opera- 
tor, one of them diverging. Nevertheless, the boundary 
approximation could be carried through, and even fully 
analytically. 

The analysis has two important merits. First, the 
curvature correction to the potential jump 0+ — (f)~ 
across the interface might remove non-regular solutions 
38] found in the "minimal regularized boundary model" 
33 . I33L l36l . I37I I . Second, it provides a very practical re- 
lation between the front curvature and the electric field 
behind the front. 

The derived boundary approximation accounts for 
propagating streamers in fairly homogeneous fields and 
far from electrodes. On the other hand, ionization fronts 
driven by monopole fields (e.g., close to sharp needle 
electrodes) or during branching undergo fast temporal 
changes of the electric field £'+ on the ionization front. 
They cannot be approximated by uniformly translating 
fronts anymore; and time dependent corrections need to 
be taken into account. In a forthcoming paper, we will 
explore how far charge conservation principles will carry 
us towards a full derivation of a moving boundary ap- 
proximation for the minimal streamer model. 
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APPENDIX A: SOLVING THE ZERO MODES 
EQUATIONS 

In this appendix, we construct the exact expression 
for the zero modes. The system of ODE we need to solve 
for the zero modes is, according to Eqs. (|63|) - ([66)) the 
following 

(/O - (7o)ljJa + Mp -tpE- {Vo + Eo)dQljJa = 0, (Al) 

o'oV'ct + V'-B - vadi^ipp = 0, (A2) 
cro/o('0(T + '4'p) - {dc,(ya)i^a + c^ = dQ%l{i^) 

Summing Eq. (UTjl with Eq. (U2|) leads to 

Eodi^tlJa = -vodcit/Ja + '4>p) + /o(V'rT + V'p)- (A4) 

We will derive below an equation for il)„ + ipp- Once we 
get it, the complete system is solved with the help of 
(jA4p and (IA3p . For this purpose, we take the derivative 
of Eq. (jA2p and replace d^4'E by its expression given by 
(jA3| to get 

crad^ipa + o-o/o(V'<T + Ipp) + C0 - vod'^tpp = (A5) 

Taking the derivative of (|A4[) and using the fact that (see 
Eqs. dSl) and (pl) l 



dcfo = -fodcEo and d^Eo 



vo 



(A6) 



we obtain 



t;o9|(V^. + ^p) - ^^^{^a + i>p) + /o5c(^a + V'p) 



Vo 



va 



dc^P^ - Eodl^„ = 0. 



(A7) 



We finally multiply (jA5|) by Eq/vq and sum it with (|A7p 
to get the equation for tpa- + ipp 

~{vo + Eo)d^^{t/jaH^p) + fodc{tlJaH^p) + ^ = 0. (A8) 

This is formally a first order ODE and can be solved 
exactly. Once 4'o- + "Pp is known, tp^ is computed with 
[4]) and then tpE is obtained from (|A3|) . 



1. Expression of il^cr + ipp 
To solve (|A8|) . we use the ansatz 



5c(V'a+^p)(C) - C(C)e^('^\ with 

5 C = / dx . 



(A9) 
(AlO) 
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The equation for C reads 

d^c= /^-^° .-g(c). 

wo(wo + £-0) 
After integration we obtain, 



where 



C(C) - C(-^) + dx^-——e 



-s(^) 



(All) 



(A12) 



^"(-) Eo{y)po{y) 



s{E^'{x))=l dy Z'^'^VT^ (A22) 



I'o + ^^oly) 



Again, with the change of variable (jA17p one can com- 
pute that 



However from the asymptotic analysis of the system of 
ODE for the zero mode (see (I75|) and ([75]) '). we know that 

C{-^) = 9c (V'. + ^p)(Cini ^ -^) = 0. (A13) 

The solution we search takes thus the form 



{Eg\x)) = dz^x 
Jq 



(A23) 



(V',+V'p)(C)=c,p+/ dye^^y^ dx 



C 



^'t'^o c-qM 



Consequently, the final form for the solution of Eq. (|A8|) 
is given by 

(V-. + V'p)(C) = c.p + ^ /''°*^' dx ^^^^, (A24) 

"0 Jo p[xy 



vo{vo + Eo) 



(A14) 



where 



c<xp = (^<T+^p) (-00). (A15) 

The explicit form of g(C) is computed using (jAlOp 



»(0 ^ /'" ..- ^»<^' 



vq + Eo{x)' 
^ £^o(x)Qo(a;) 
_oo vo + Eo{x) ' 



Using the change of variable 

y = £^0(3^) 



(A16) 



(A17) 



dy = da;£:o(a;)d2; = „ . . fl;/, (A18) 

Wo + £^o(a;) 



we obtain 
5(0 = 



iJo(C) 



dx 







ai\x\) _ /■'^°'''^' a{x)_ 
P{x) Jo P{x) ' 



p(x) vp(|£^o(c)I); 



P(0) 



In 



(7 



Po(C) 



(A19) 



where we use the notation p(|i?o(C)|) = p{Eo{0) — PoiC) 
in this Appendix. Substituting Eq. (|A19D in Eq. (|A14D 
we obtain 



Poiy)J-oo Vo{vq+Eq) 



With the change of variable (|A17p we obtain 



(A20) 



(^, + ^p)(C) = C.P + - J^ dx -^-p , 

(A21) 



2. Expression of ^a 

We have now to solve Eq. (|A4p . We rewrite it using 
the exact solution (|A24[) . In particular we have 



9c(V'cr + V'p 



Thus 



£^!!ii_^aci?o 

"0 Po 

^0 Po 



dci^a = ao{tpa + Ipp) 

Po 



(A25) 
(A26) 

(A27) 



Remark that it is easy to verify with the help of this re- 
lation that we can recover the asymptotic behavior (|75p . 
After an integration of (IA27P we obtain 



V'cr = V'(T(Cini) -C0 



dxao{x) 



dx 
Cini ^0(2;) 

«o Jo p[y) 



From the asymptotic analysis of (jA27[) . and as we show 
explicitly below, one can verify that only the first integral 
of this last expression diverges when (^{^1 — >■ —00. Using 
the change of variable (IA17P we finally get 



V'ct(C) = '0<T(Cini) -C0 



\Eo{0\ 



dx 



Vq- X 



"•0"P 



l^«(^)l ^ {vo-x)a{x) 



+ 



dx- 
dx- 



|Bo(Cini)| 2;p(a;) 

C0 /■'^"('^^l _^^(f o-x)a(x) 
xp{x) 



^0 J\Eo{UO\ 



Vq 



piy)' 



-29) 



In this last expression only the first integral diverges 
when Ciiii — > —00. The singularity can be explicitly ex- 
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tract thanks to the relation 



dx 



Vq - X 



+ 



|£o(C.„i)l ^Pi^) 

C - Cini 



^°^'^^' , VQ-X ^ Cini - C 

dx — -—^ H 

|So(OnOI XP[XV CT 



a 

I So (01 



da; 



Wo -x 



|£o(Ci„i)l 2;p(a;) 

C - Cini 



va~ X 



So (01 

da;- , , 

Bo (Cini) I ^^(2^) 



|So(OI 

I Bo (Cini) I 



da; 



1 



p{x) cr- 



^'o -.^ C-Cii 
a;p(a:) cr^ 



-, (A30) 



where we used Eq. ([22]) ■ The integral of this last relation 
is always finite. Expression (jA29p can then be written as 
in the main text 



V'.(C) = ^Cini + a.(C)+&.(C)e^^" 
a 



where 

^^(0 = ^^(Cini) -C0 



Bo(OI 



da; 







p{x) a- 



(A31) 



xp{x) 



_ 2^, /■l^°(^)l (^;o - x)a{x) 

a Jq xp(x) 

vo Jo xp[x) Jo p[yy 

The expression 6o-(C) is not useful for our analyze since 
the contribution of this term in Eq. (jA3ip can be made 
as small as we want since jCinil can be as large as we want 
(but finite). 



3. Expression of ipE 



Equation (jASp can be written as 

C^CV'-E = cro/o(V''cr + i'p) - (9cf^o)V'(T + C0 

= cro/o(V''cr + V-'p) + o-od(ipcr - di^iaoipa) + c^ 

= 0-o/o(V^ff + V^p) - 0"o^ CToQ:o(V^<t + V'p) 

Po 



C^C('^O'0cr) +C0 



(A33) 



where we used the expression (jA27|) for dQi/ja- Using the 
relation ([45]) for /q, this last expression reduces to 

dc'ipE = — F-^(V'<T + V'p) ^-4r - dc{(jQtlJa) + C0, 

-C^O VO + -C/0 

(A34) 
where we used also the relation (1201) between ctq and pQ. 
It is also easy to see that in the limit C ^ — oo we re- 
cover the asymptotic form (j77p . The expression of ipE is 



obtained after an integration of (|A34p 

^£;(C) = V'BlCini) 



ax c„ 



Eo 

«o Jo p{y? 

+ Cro(Cini)l/'cr(Cini) - C^Wq 



-^0(0^-^(0 

'^ da; 



Ci. 



vo + Eo 

+ C0(C- Cini) (A35) 

Using the change of variable (jA17[) we get 



V'b(C) 



V'BlCini) - Vo I dx 5— [Cap 



+ ^ I dy 



So (Cini) 

(«o + y) 



1^0 Jo " p{yY 

O-Q ( Cini )V'o- (Cini) - C^Uq 



-^o(C)^^(C) 

^"(^ da; 

S0(C>ni) ^Pi^) 



+ C^C-Cini) (A36) 

The two last terms diverge when Cini — > — oo but their 
sum is finite. With the help of the definition (|22p we can 
write this sum as a finite quantity. Indeed, we know that 

f-So(C) 

C - Ci. 



vo+ X 



dx- . , 
So(C.„0 ^P^^> 



(A37) 



and we substitute it in (jA36| to get 

r\Eo{C.)\ Q,(-^^ 

V's(C) = V'B(Cini) + CapVo / 2~ 

- C0 / c^a;^^ / dy — -^-^ cro(C)'0<T(C) 



CTO ( Cini )V'o- (Cini) - C0 



piyf 

^«(«)l dx 



(A38) 



/|So(C.nO|/'(^)' 

The only part of this relation that diverge when Cini -^ 
— cx) is contained in the expression of il^aiC)- We can still 
simplify this last relation by using the following formulas 

i|So(0„OI -^ ~ ^"^"^J|-0(Cini)| 

= a(|i?o(C)l)-«(|i?o(Cini)|).(A39) 

An important point to emphasize is that this calculation 
is the only one so far where the explicit form of the func- 
tion a{x) = e^^'^ was used to simplify the expression. 
The second relation is obtained after an integration by 
parts 

r"'da;^rd,-(^°+^)- 



p{vf 



r\Eo{C)\ 
{a{\Eo{C)\) ~ a{\Eo{Cini)\)) dx 

"'lSo(Cini) 

^»('^)l , a(x)(vo-x) 
da;- 



|So(OnOI 



p{x) 



Vo - X 

p{xY 

(A40) 
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Finally we get 

■0£;(C) = c^ini +Cap«o(ao(C) -ao(Cini)) 



C0(ao(C) -ao(Cini)) 
dx 



So (01 



dx 



vo- X 



C0 



C0 



^"(^^1 , a(a;)(wo-a;) 



|iJo(C)l rf^ 



|So(Ci„i) 



p{x)' 



'To(C)V'.(C) 

(A41) 



with 



CC.n. = ^EiCini) + MCin O^aiCini) (A42) 

Expression (|A4ip can then be written as in the main text 



a 



where 



\Ea(Q\ 



(A43) 



vo~x 



aEiC) = c^i,,i + CapVoaoiC) - C4,ao{C) / dx 

Jo P\^ 

r\Eo(Q\ a{x){vo-x) , ^, .,^,,, 

+ c^j dx ^ ^y^^ ^ -ao(C)ai')(C) 



C0 



^"('^^l dx 



(A44) 



/o P(2^)' 

where we used (|A31I) and where 

cCini = V'sCCini) + o--i/;„(Ci„i). (A45) 

As for (jA3ip . the expression 6b (C) is not useful for our 
analyze since the contribution of this term in Eq. (jA43p 
can be made as small as we want. 



4. Expression of A 

Since we know that Eq, ctq and po decay exponentially 
for ^ — + — cx), see (^5)) . and since the zero modes diverge 
at most linearly in this limit, we can already replace C,c 
by — oo to compute A and B. 

We know the exact expression of the zero modes and we 
start to compute the relevant quantities for the relation 
(j62p which will give us the boundary conditions we search 
for. We use in this section the notation 



(t/-, +7/v)(C) - C„p + C^h{EQ{0), 



(A46) 



h{E,{C)) = - r^^\jj^^l±A, (A47) 

Wo Jo pixr 

From the definition ([57|l of A and Eq. (JA46| , we have 



A 



dC bPadi^a-Q + (ccrp + c^h - ijjc) 5f po] 



-CT C„p+ d(lpaid((7o- d(;po) 



+ c^ dChdi^po 



(A48) 



Using p4|) we can compute the derivation of the ion den- 
sity of the planar front solution 



a a fvQ + Eo 

V ^0 
_ E'pgg vo + Eq 
vl Vo 



Oqcto. 



We substitute it in Eq. (|A48p to obtain 



A 



-a c. 



ap 



dCtP, 



Eocri 



dC 



+ c^ dChdt^po 



(A49) 



^paEQ^Qao 
Vo 

(A50) 



The two last integral can be simplified. Indeed, using an 
integration by parts we have 

/ dC,tp^EodQ(To = [ilj„Eocrof_^ 

f' — oo 



— oo 




dC (To {EodQi>a + ■Ipad(Eo) 



(A51) 



where we have still to replace d^ipa and d(^Eo by their ex- 
pression (see (|X27| and ^). The last integral of (|M0)) 
can also be simplified by an integration by parts 



dChdQpo 



dC- 



Eq 

Vo 



(A52) 



Substituting Eqs. (|A5ip and (|A52p in Eq. (|A50p . making 
the change of variable (jA17p and using the definition 
of (7^, we get 



A = cro(O)V'„(O) + C0 







Ctjj I dxa{x)h{—x) ^ — - i dx 
/o ^0 Jo 



dx 
p{x) 

C0 



(A53) 

vo-x 
p{x) 



Notice that the only term that diverge when C,^-^^ — > — cx) 
is contained into ^a{Q). Now this formula can be simpli- 
fied further since actually the two last integrals cancel. 
Indeed, with an integration by parts we get 



C0 / dxa{x)h{-~x) 

10 

E^ 



C0 
Vo 



H{x) f dy 
Jo 



vo-y 



lE+ 



p{y? 



+ ^/ dx\^{~pix)+H{E+)), 

Vo Jo pyxy 

where we used the notation 



(A54) 



H{x) = / dxa{x) = - / dya{y)+H{E+), 

= ~p{x) + H{E+). (A55) 

It is easy to verify that H{0) = 0, consequently we have 

Cd, / dx a{x)h(—x) = / dx — r^— ■ (A56) 

Jo vq Jq p{x) 
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The final expression of A is then rather simple and reads 



A = cro(O)V'a(O) + C0 



dx 



(A57) 



/o Pix) ' 
This expression can be written as in the main text 



where 



A^^^^Cin^ + aA + bAe'<'-, (A58) 



a A = cro(0)ao.(0) + c^ 



dx 
P{x) 



(A59) 



5. Expression of ^ + V'-b(O) 



We just need to sum the expression (|A57p of A with 
the expression (jA41[) of ■i/'s evaluated at C = 0. We get 



^ + V^£;(0) = c^i„ + c^pWo(a(^+) - ao(Cini)) 
- C0(a(£'+) - ao(Cini)) / dx- 

J\Eo{UO\ 



p{xY 



C4, 



E+ 



dx 



|£0(C.ni)| 



a{x){vo -x) 
p{xY 



(A60) 



with Ci^i^j defined by (|A42p . This expression can be writ- 
ten as in the main text 



A + Vb(0) = oa + 0^(0) + 6B(0)e^';'"', (A61) 



where 



dA + aB(0) = c^i„i + CcrpVQa(E^) 
C0a(i?+) / dx ^^ 



vq-x 

p{xy 



a{x){vo — x) 
P{x) 



+ Cci, I dx ^^ 



(A62) 



and where c^^^^^ is defined by (jA45[) . 



Replacing ipE by its exact form (|A41[) we get 

B = -j dx — ——[cQ,^_^.^+CapVo{a{x)-a{\Cini)\) 

JQ P(X) 

^° -pix)^^ {E-\-x)) 



Vo- X 

C4>{a{x) -a(ICini)l) 



vp-y 



dy-^,^.2 



, , a.{y){vo -y) 

Now we need to compute 

,E+ 

Vo-X 



Bo(C.„i)l ^12^ 



c,p I dy 



'0 J|Bo(OnOI XP{X)^ 

- c., /"' dy r dx ^^'^^ff^ (A65) 

+ ^ /''" dy r dx ("° -")"(") /' dz^^ 

^'0 7o JiBoCOnOI ^/'(^) -'o P[zy 

The repeated integrals can be computed by parts to get 

rE+ 

dxipa {Ei^'^{-x)) = '(/'<T(Cmi)t'0 


c^-wo / dx — —-:r + C0 / dx 

-'|Bo(Ci„i)l a;p(x)^ 7o P[xy 

/■^ (i;o - a;)Q;(x) 

CapVo / ax ;r^; 







(uo ~ x)a{x) 
p{x) 



(A66) 



C0 



^^(.o-x)a(x) f\yVo-y 



|£0(Cini)| 



xp{x) 



pivY 



C0 /"^ [vq " a::)a(a::) r ^^VQ-y 

va Jo P{x) Jo 



''-piyy 



For the first integral in Eq. (|A66p . we use Eq. (|A30p in 
order to extract the part which diverges as Cini ~> —00. 
Now we substitute Eq. (|A66p in Eq. (|A64p . and in the 
final result, in order to eliminate a maximum of double 
integral, we use the following relation 



6. Expression of B 



Using the change of variables (jA17p in the definition 
of B we get 



B 



E+ _ 

dx^T-^i'E {Eo\-x)) 

pyx) 



(A63) 



P{x) J\Eo{UO\ P^vY 



dx- 



dx 



P{x)' 
vq-x 
P{x) 



|So(C.„OI 



Po(Cini) JO 

dy 

p{y) 



dx 



p{x) 
(A67) 
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The final expression for B is tlien given by 

CT Jo P[X) 



Tfiis expression can be written as in the main text 



c^vl I dx 



2 

B = ^^Cini + aB+ bee^'^" 



(A69) 



1 1 



p{x) (j- 



E+ 



xp{x) 

E+ 



[vq - x)a{x) 



where 



/ 7 U -^ 2 1 7 ^ 

+ c^vo / dx - CapVQ / dx- 

Jo p{xy ■/|Bo(Ci„i)i ^py^i 

[^ 7 («o - x)a{x) r (wo - y) 
+ C0Wo / dx -J— / dy — --2— 



as = WoV'o-(Cini) -c^ini / da; 



woj-a; 
p{x) 



dx 



Vq - X 



Po(Cini) Jo Pi^) 

E+ 



f 

c^Vq / dx 
Jo 

rE+ 



1 1 



p{x) a- 



; 7 wo-x^ /!/• i\ / J ^0 -a; + c^wo / dx - c^pV^ dx 

c^ dx --^ h c^pWoa( Cini ) / dx — —— Jq p{xY Jq xp{x) 

Jo P{xY Jo p{x) 



xp(x) 

rE+ 



pi^y 



n. l^ I 7 (^0 - x) f {vq - y) 

C0a(Cini ) / dx —— / dy — -—^ 







/o P{x) 



dy 



B0(Cini)| 



(^'o - V) 

p{yf 



I J (^0 - x)a{x) r {vq - y) 
"'^^""/o '^ Mx) L'^'l^ 



+ ^ / dx ^^ 
CT Jo p{x) 



C0 / dx- 



p(x) 



,(A70) 



I J ("0 - a;)a(a;) T (uo - y) 



(A68) with a^^j^j defined by Eq. (g45|. 
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